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Abstract. Fhe interferometric observations of dwarf stars in the solar neighbourhood, combined with HlPPARCOS 
parallaxes provide very precise values of their linear diameters. In this paper, we report the direct measurement 
of the angular diameter of the bright star Sirius A with the VINCI/VLTI instrument. We obtain a uniform disk 
angular diameter of #ud = 5.936 ± 0.016 mas in the K band and a limb darkened value of #ld = 6.039 ± 0.019 mas. 
In combination with the HlPPARCOS parallax of 379.22 ± 1.58 mas, this translates into a linear diameter of 
1.711 ± 0.013 D0. Using the VINCI/VLFI interferometric diameter and the published properties of Sirius A, 
we derive internal structure models corresponding to ages between 200 and 250 ± 12 Myr. Fhis range is defined 
mainly by the hypothesis on the mass of the star, the overshoot and the metallicity. 
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1. Introduction 

The brightest star in the sky, Sirius A (HR 2491, HD 
48915), is an A1V type dwarf. It is also a member of 
the fifth nearest binary system. The existence of a mas- 
sive companion, Sirius B, was predicted by Bessel l|1844f) 
based on the apparent motion of Sirius A. It was later 
discovered visually by A. G. Clark in 1862, and is still to- 
day among the most massive white dwarfs (WD) known. 
Chandrasekhar (4935) published the first model to explain 
its characteristics. Sirius has been reported several times 
to be a red star in the ancient times, whereas it is now 
almost white. The reason behind this color change is still 
not understood, though several scenarios have been pro- 
posed in the last decades (Schlosser & Bergmann [1985 
Bonnet-Bideau & Grv ITMTl Whittet H^jl . During the 
past century, many efforts were directed towards compil- 
ing the properties of both stars like their masses or chem- 
ical abundances. Sirius A exhibits enhanced metal abun- 
dances in its photosphere and has therefore been classi- 
fied as Am by Abt 1^)79). As the orbital period is rela- 
tively short (50 years), two revolutions of the Sirius A-B 
binary have already been measured, allowing Gatewood & 
Gatewood (|1978[) to derive precisely the masses of the two 
components (Sect. l2~T)l . Moreover, Sirius A has been the 
first main sequence star to have its diameter measured in- 
terferometrically by Hanbury Brown ct al. ( 19671 lHj74a[) . 
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The chief object in our study is to model the internal 
structure of Sirius A in order to reproduce its macroscopic 
characteristics. In this process, we use the angular diame- 
ter of this star from our new VINCI/VLTI measurements 
and the HlPPARCOS parallax to derive its linear diame- 
ter. From our models, we compute the age and the initial 
helium content of Sirius A. We summarize in Sect. [21 the 
fundamental parameters of this star, while Sect. |3| to 
are dedicated to the VINCI data processing and analy- 
sis. This section includes in particular the fit of the limb 
darkened disk model to the measured visibilities. Finally, 
we discuss in Sect. the macroscopic constraints that we 
apply to our CESAM numerical models (Morel IT99*7|) of 
the internal structure of Sirius A, and we predict the as- 
teroseismic large frequency spacing of this star. 

2. Fundamental stellar parameters 

2.1. Masses 

Sirius (A1V + WD) is a visual binary system with a 
50 years period that has been studied extensively, and 
its orbit is known with high precision. From photomet- 
ric astrometry covering a 60 year period, Gatewood & 
Gatewood (|1978f) determined the orbital photocentric 
semi-major axis of Sirius A with respect to the system 
barycenter. This determination relied upon the orbital 
parameters of the relative visual orbit of Sirius B given 
by van den Bos (1960). These authors derived a mass 
for Sirius A of 2.143 ± 0.056 M Q , using a parallax of 
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7r = 377.7 ± 3.31 mas. This precision of tt has been re- 
cently improved thanks to the Hipparcos satellite mis- 
sion (Perryman et al. ITW7|> at tt = 379.22 ± 1.58 mas. 
We note that this value is consistent with the paral- 
lax of Van Altena et al. {1995), tt = 381.6 ± 2.2mas. 
Consequently, the sum of the masses has to be dimin- 
ished by 1.2%, giving for Sirius A a lower mass of 2.12 
± 0.06 Mq. Independently of binary orbital elements, 
Holberg et al. ((1998 ) have re-determined the mass of 
Sirius B from precise measurements of its surface grav- 
ity and radius. Using the Hipparcos parallax they found 
a mass of Mb = 1.034 ± 0.026 M Q in agreement with the 
revisited sum of the masses. Also based on the Hipparcos 
parallax, Provencal et al. ( 1998) have derived a value of 
M B = 1.000 ± O.O16M , in agreement with Holberg et 
al. lfTM5|l at a 1.1 a level. 

2.2. Age 

The age of the Sirius system has been poorly debated in 
the literature. Based on Wood's (1995) sequences of cool- 
ing age of WD (carbon-oxygen model with thick H layer) , 
Sirius B is found to be a young WD of about 160 (Holberg 
et al. H998[) to 210 ± 20Myr (Holberg, private communi- 
cation). The revisited mass of Sirius B (1.034 M ) corre- 
sponds to a progenitor mass of « 7M Q , i.e. the approxi- 
mate mass of a B5 main sequence star, from the empirical 
initial to final mass relation of Weidemann (2002J. Such a 
massive star evolves to a WD in 40 ± 5 Myr. Consequently, 
assuming a simultaneous origin of both A & B compo- 
nents, the Sirius system is expected to be between 200 
and 250 Myr old. 

2.3. Spectral properties 

Sirius shows deep metallic spectral lines, and we propose 
that the corresponding enhancement in metal abundance 
is a consequence of the strong radiative forces at work 
in the interior of this star. This hypothesis is supported 
by the very low rotational velocity (Royer et al. I2002f) 
vs'mi w 16 km s , that favors the action of these forces. 
In their comparison with the star Vega, Qiu et al. (2001) 
concluded that Sirius A appears to be enhanced in abun- 
dances by about 1.0 dex on average while C, O, Sc and 
Ca are underabundant. In consequence, for such a slow 
rotator and intermediate mass, it is very hazardous to es- 
timate the mean metal content in the interior of the star 
from its surface metal abundance. 

Therefore, we propose to focus only on the abundances 
of CNO, that arc believed to be little modified by radiative 
forces (Richer et al. EDOOi), to estimate the value of (-§■) 
at present day. We adopt an average abundance at present 
day of (-§)_, = 0.012 ± 0.002 in our models, that includes 
the microscopic diffusion of the elements. (-§-) ; is very 
important for the internal stellar opacity of the computed 
model and is a key parameter with respect to the age of the 
star (see Section EJ). It results, when using these elements 



as tracers, that the average metal abundance of Sirius A 
is in fact about half that of the Sun. While appearing 
to be low, the adopted Z s abundance is comparable with 
Vega's observed value. Vega is a fast rotator and therefore 
should not be affected by radiative accelerations masking 
its true metal abundance. This corresponds to an initial 
value at zero age of = 0.0165 ± 0.002. In order to 

validate the previous assumption that the metallicity of 
outerlayers does not represent the metallicity of the whole 
star, we also present, in Sect. El a model of Sirius A with 
the observed metallicity (y) s = 0.0329, corresponding to 
the observed mean surface abundances. 

2.4. Photometry 

The photometric properties of Sirius are summarized in 
Tabled The magnitude and the parallax give an intrinsic 
luminosity of L/L Q = 25.4 ± 1.3. Following Lemke ifH^jl . 
we adopt an effective temperature of T c ff = 9900 ± 200 K 
and a logarithm of the spectroscopic surface gravity of 
\ogg = 4.3 ± 0.1. 

3. Interferometric observations and data 
processing 

3.1. Instrumental setup 

The European Southern Observatory's Very Large 
Telescope Interferometer (Glindemann et al. I2000f) is op- 
erated on top of the Cerro Paranal, in Northern Chile 
since March 2001. For the observations reported in this 
paper, the light coming from two test siderostats (0.35 m 
aperture) was recombined coherently in VINCI, the VLT 
INterferometer Commissioning Instrument (Kervella et al. 
I3UUU|) . We used a regular K band filter (A = 2.0 — 2.4,um) 
but VINCI can also be operated in the H band (A = 
1.4 — 1.8/rai) using an integrated optics beam combiner 
(Berger et al. EMTl Kervella et al. I2003ajl . Two VLTI 
baselines were used for this program: D1-B3 and E0-G1, 
repsctively 24 and 66m in ground length. 

3.2. Data processing 

We used a modified version of the standard VINCI data 
reduction pipeline (Kervella et al. 2003c), whose gen- 
eral principle is based on the original FLUOR algorithm 
(Coude du Foresto et al. 1997). The two calibrated output 
interferograms are subtracted to remove residual photo- 
metric fluctuations. Instead of the classical Fourier analy- 
sis, we implemented a time-frequency analysis (Segransan 
et al. 1199 9;) based on the continuous wavelet transform 
(Fargc 1992). In this approach, the projection of the signal 
is not done onto a sine wave (Fourier transform), but on a 
function, i.e. the wavelet, that is localized in both time and 
frequency. We used as a basis the Morlet wavelet, a gaus- 
sian envelope multiplied by a sine wave. With the proper 
choice of the number of oscillations inside the gaussian 
envelope, this wavelet closely matches a VINCI interfero- 
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gram. It is therefore very efficient at localizing the signal 
in both time and frequency. 

The differential piston corrupts the amplitude and the 
shape of the fringe peak in the wavelets power spectrum. 
A selection based on the shape of fringe peak in the time- 
frequency domain is used to remove "pistonned" and false 
detection interferograms. Squared coherence factors /j, 2 are 
then derived by integrating the wavelet power spectral 
density (PSD) of the interferograms at the position and 
frequency of the fringes. The residual photon and detector 
noise are removed by making a least squares fit of the PSD 
at high and low frequency. The resulting measurement sta- 
bility is satisfactory: on a good series of 500 interferograms 
of Sirius (10 minutes) the standard deviation is typically 
2.0%, and the final bootstrapped statistical error on the 
average fj, is only 0.16%. The /x 2 values are converted into 
calibrated visibilities V 2 through the observation of cali- 
brator stars (Sect.UJ. 

4. Calibration 

The calibration of Sirius visibilities was achieved using 
well-known calibrator stars that have been selected in 
the Cohen et al. ( 1999 ) catalogue. The characteristics of 
the selected calibrators are listed in Table ^ The limb- 
darkened disk (LD) angular diameter of these stars was 
converted into a uniform disk value using linear coeffi- 
cients taken from Claret et al. I|1995[) . As demonstrated 
by Borde et al. I|2002[) . the star diameters in the Cohen 
et al. (1999) list have been measured very homogeneously 
to a relative precision of approximately 1% and agree well 
with other angular diameter estimation methods. The in- 
tcrfcromctric efficiency is computed from the coherence 
factors obtained on these calibrators, taking into account 
the bandwidth smearing effect (Sect. 15. 2|) and a uniform 
disk angular diameter model. This process yields the cali- 
brated squared visibilities that are used for the model fit. 
In Table El two error bars are given for each V 2 value: 

— one statistical error bar, computed from the dispersion 
of the visibility values obtained during the observation, 

— one systematic error bar defined by the uncertainty on 
the knowledge of the calibrator stars angular sizes. 

While the statistical error can be diminished by repeat- 
edly observing the target, the systematic error cannot be 
reduced by averaging measurements obtained using the 
same calibrator. This is taken into account in our model 
fitting by checking a posteriori that the resulting uncer- 
tainty of the model visibility is larger than the systematic 
errors of each measured visibility value. This conservative 
approach ensures that we are not underestimating the fi- 
nal error bars on the angular diameter of Sirius. 

5. Angular diameter computation 

5.1. Intensity profile 

From the visibility curve shape before the first minimum, 
it is almost impossible to distinguish between a uniform 
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Fig. 1. Overview of the VINCI/VLTI visibility measure- 
ments obtained on Sirius. The solid line is the best fit vis- 
ibility model, that takes into account the limb darkening 
and the bandwidth smearing. 
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Fig. 2. Visibility measurements obtained on Sirius on the 
E0-G1 baseline. The thin line is the best fit model, and 
the dashed lines mark the limits of the la error domain. 
For clarity, the plotted error bars are the quadratic sum 
of the statistical and systematic errors. 

disk (UD) and limb darkened (LD) model. Therefore, it 
is necessary to use a model of the stellar disk limb dark- 
ening to deduce the photospheric angular size of the star, 
from the observed visibility values. The intensity profiles 
that we used were computed by Claret l|2000|l . based on 
model atmospheres by Kurucz l|1992ll . We chose the four 
parameters /(/j,)//(1) law of this author, where /i = cost9 
is the cosine of the azimuth of a surface element of the 
star. They are accurate approximations of the numerical 
results from the ATLAS modeling code. 

To compare the different types of approximated limb 
darkening laws, we chose the following approximate pa- 
rameters for Sirius: T cS = 10000 K (for Claret et al.[l995 
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Table 1. Relevant parameters of aCMa and its calibrators from the literature. 



Name 


aCMa 


(9Cen 


8 Lep 


aCrt 


310ri 


HR 4050 


HD number 


HD 48915 


HD 123139 


HD 39364 


HD 95272 


HD 36167 


HD 89388 


mv 


-1.47 


2.06 


3.80 


4.07 


4.71 


3.38 


rriK 


-1.31 


-0.10 


1.31 


1.62 


0.90 


0.60 


Sp. Type 


A1V 


KOIIIb 


G8III/IV 


K1III 


K5III 


K3IIa 


Tea (K) a 


9900 


4980 


4580 


4650 


3930 


4500 


[Fe/H] f 


0.5 


0.0 


-0.75 


-0.2 


-0.3 


0.5 


logs' 


4.3 


2.75 


2.95 


2.8 


1.6 


1.6 


7r (mas) b 


379.22 ± 1.58 


53.5 ± 0.8 


29.1 ± 0.6 


18.7 ± 1.0 


7.2 ± 0.8 


4.4 ± 0.5 


#ld (mas) c 


6.039 ± 0.019 s 


5.434 ± 0.020 e 


2.63 ± 0.04 


2.28 ± 0.03 


3.66 ± 0.06 


5.23 ± 0.06 


#ud (mas) d 


5.936 ± 0.016 g 


5.305 ± 0.020 c 


2.57 ± 0.04 


2.22 ± 0.03 


3.55 ± 0.06 


5.09 ± 0.06 



a From Cohen et al. (1199211 for a CMa, and from Cayrel de Strobel et al. (119971 1200H for other stars. 

b Parallaxes from the HlPPARCOS catalogue (Perryman et al. 119971 . 

c Catalogue values from Cohen et al. (19991 for other stars. 

d Linear limb darkening coefficients factors from Claret et al. (1995). 

° OCen angular diameters from Kervella et al. 12003bt . 

f Qiu et al. for a CMa, and Cayrel at al. 1119971 for other stars. 

g Sirius angular diameters from this work (see text for details). 



Table 2. Sirius squared visibilities. 



JD 




Stations 


Baseline (m) 


Az. b 




V 2 


± Stat. ± Syst. 


Calibrators 3 




2452682 


.71653 


D1-B3 


17.758 


81.68 


0. 


8602 


± 


0.0206 


± 


0.0007 


8 Lep, 


aCrt 




2452682 


.71251 


D1-B3 


18.139 


81.45 


0. 


8742 


± 


0.0210 


± 


0.0007 


8 Lep, 


aCrt 




2452679 


.71932 


D1-B3 


18.270 


81.36 


0. 


8769 


± 


0.0313 


± 


0.0023 


8 Lep, 


31 Ori, 


HR4050 


2452682 


.70813 


D1-B3 


18.544 


81.18 


0. 


8668 


± 


0.0208 


± 


0.0007 


8 Lep, 


aCrt 




2452679 


.71478 


D1-B3 


18.684 


81.08 


0. 


8748 


± 


0.0309 


± 


0.0023 


8 Lep, 


31 Ori, 


HR4050 


2452679 


.71085 


D1-B3 


19.032 


80.84 


0. 


8627 


± 


0.0300 


± 


0.0022 


8 Lep, 


31 Ori, 


HR4050 


2452679 


.66856 


D1-B3 


22.063 


78.01 


0. 


8092 


± 


0.0280 


± 


0.0021 


8 Lep, 


31 Ori, 


HR4050 


2452679 


.66415 


D1-B3 


22.301 


77.68 


0. 


8110 


± 


0.0281 


± 


0.0021 


8 Lep, 


31 Ori, 


HR4050 


2452679 


.65967 


D1-B3 


22.527 


77.35 


0. 


8048 


± 


0.0278 


± 


0.0021 


8 Lep, 


31 Ori, 


HR4050 


2452361 


.62629 


E0-G1 


60.439 


1.29 


0. 


1599 


± 


0.0074 


± 


0.0045 


0Cen 






2452361 


.64154 


E0-G1 


60.574 


5.16 


0. 


1608 


± 


0.0062 


± 


0.0045 


OCen 






2452361 


.57907 


E0-G1 


61.038 


169.41 


0. 


1472 


± 


0.0085 


± 


0.0041 


OCen 






2452358 


.55861 


E0-G1 


62.027 


162.72 


0. 


1433 


± 


0.0062 


± 


0.0042 


6Cen 






2452340 


.60484 


E0-G1 


62.148 


162.07 


0. 


1311 


± 


0.0059 


± 


0.0014 


OCen 






2452340 


.57040 


E0-G1 


63.651 


155.00 


0. 


1243 


± 


0.0063 


± 


0.0014 


OCen 






2452340 


.56629 


E0-G1 


63.831 


154.23 


0. 


1226 


± 


0.0054 


± 


0.0013 


6Cen 






2452340 


.54577 


E0-G1 


64.678 


150.65 


0. 


1186 


± 


0.0051 


± 


0.0013 


OCen 






2452340 


.54210 


E0-G1 


64.816 


150.06 


0. 


1152 


± 


0.0052 


± 


0.0013 


OCen 







a All calibrators for Sirius were chosen in the Cohen et al. l|T|599} catalogue. 
The azimuth is counted clockwise from North, in degrees (E = 90 deg). 



and Claret I2UUU|) or T cff = 9800 K (for Claret EMS) , 
log 5 = 4.0, [Fe/H] = +0.5, V tmh = km.s" 1 . These were 
chosen as close as possible to the Cohen et al. (1992) val- 
ues. As noted by Claret l|2000J) . the sensitivity of the limb 
darkening of a T c tf — 10000 K star to the metallicity is 
negligible, so the true metal content of Sirius is not crit- 
ical. Different I(/i)//(l) square root and four parameters 
laws for Sirius, from Claret et al. (1995) and Claret (1998 
2000| are plotted on Fig. EH Though no formal error bars 
are provided for these limb darkening models, the small 
differences between the four curves demonstrates the good 
internal consistency of the different types of laws. 



The limb darkening is directly measurable by inter- 
ferometry beyond the first minimum of the visibility 
function, as demonstrated by several authors on giant 
stars (Quirrenbach et al. ] 1996 , Wittkowski et al. I20()l|l . 
Unfortunately, for stars of the angular size of Sirius ob- 
served in the K band, this requires baseline of 80 to 120 
meters that were not available for the measurements re- 
ported here. It is intended in the near future to measure 
directly the LD of a number of nearby stars, using the 
VINCI and AMBER (Petrov et al. l2TOr]| instruments on 
the long baselines of the VLTI. With lengths of up to 
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95% 



. 85% 



70% 



sqrt 2000 
- 4par 2000 
sqrt 1998 
sqrt 1 995 



Fig. 3. Published LD laws for Sirius in the K band, from 
Claret et al. l(TM5jl . Claret (|P39g|> and Claret pUUUjl . The 
square root law is already a good approximation to the 
Kurucz atmosphere models. The solid line shows results 
with the most recent four parameters law introduced by 
Claret ( 2000). This last version was chosen for the Sirius 
LD fit of this paper. 

202m, they will allow the exploration of the secondary 
and higher order lobes of the visibility function. 

5.2. Visibility model 

The VINCI instrument bandpass corresponds to the K 
band filter (2-2 .4 /im). An important effect of this rela- 
tively large spectral bandwidth is that several spatial fre- 
quencies are simultaneously observed by the interferome- 
ter. This effect is known as the bandwidth smearing. It is 
usually negligible for V 2 > 40%, but this is not the case 
when visibilities come closer to the first minimum of the 
visibility function. In the case of Sirius A observed with 
the E0-G1 baseline, V 2 is already low ( 13%) and we have 
to consider it. 

To account for the bandwidth smearing, the model 
visibility is computed for regularly spaced wavenumber 
spectral bins over the K band, and then integrated to 
obtain the model visibility. We assume in this paper a 
limb darkening model following the four parameters law 
of Claret pWill : 



/(/*)//(!) = 1 -£>*(! 



(1) 



fe=i 



The a^ coefficient are tabulated by this author for a wide 
range of stellar parameters (T^ff j l°g<?v) an( i photometric 
bands (U to K). For Sirius (V T = km/s, T cS = 10000 K, 
log # = 4.5, [M/H] =0)wefmdai = 0.5318, a 2 = -0.6921, 
03 = 0.6698, and a^ = —0.2431. Fig. shows the corre- 
sponding intensity profile In this paper, we as- 
sume that the limb darkening law does not change over 
the K band, which is reasonable for a hot and compact 



stellar atmosphere, but this computation can easily be ex- 
tended to a wavenumber dependent a) intensity pro- 
file. Following Davis et al. (2000 ), using a Hankel integral, 
we can derive the visibility law V(B, #ld, c) from the in- 
tensity profile: 



1 f 1 

- J lit 1 ) Jo (ttB<t8l D \J\ - fjp^fi dn 



V =A 



where a is the wavenumber: 
a = 1/A 

and A is a normalization factor: 

A = / I(fj,)/j,diJ, 
Jo 



(2) 



(3) 



(4) 



The integral of the binned squared visibilities is com- 
puted numerically over the K band and gives the model 
V 2 for the projected baseline B and the angular diameter 
#ld through the relation: 



V 2 (6 hD ,B)= ( [V(B,9 hu ,cr) T(a)} 2 da 
Jk 



(5) 



where T(a) is the normalized instrumental transmission 
defined so that 



/ T(cr)da 



1 



(6) 



We computed a model of T(cr) by taking into account 
the instrumental transmission of VINCI and the VLTI. 
It was first estimated by considering all known factors 
(filter, fibers, atmospheric transmission,...) and then cali- 
brated on sky based on several observations of bright stars 
with the 8 meter UTs (see Kervella et al. 2003b for more 
details). This gives, for Sirius, an average measurement 
wavelength of 2.176 fim. 

The V 2 (6i i d,B) model is adjusted numerically to the 
observed (B, V 2 ) data using a classical x 2 minimization 
process to derive #ld- The reduced \ 2 of the global fit 
is 0.30, characteristic of a satisfactory correspondence be- 
tween the model and the measured V 2 values. Fig. ^ and 
[5] show the position of the VINCI points on the squared 
visibility curve of the chosen model. 

5.3. Angular and linear diameters 

The \ 2 minimization gives a limb darkened angular diam- 
eter #ld = 6.039 ± 0.019 mas for Sirius A, while a simpler 
uniform disk model yields #tjd = 5.936 ± 0.016 mas. The 
conversion formula between the LD angular diameter and 
the linear diameter D (in solar unit) is: 



D = 



2p tan(f? /2) 



(7) 



Assuming a solar angular radius of Oq/2 = 959.64 ± 0.02" 
(Chollct & Sinceac 1999 ) and the HiPPARCOS parallax p = 
379.22 ±1.58 mas, we obtain Dshius = 1.711 ± 0.013 D©. 
It should be noted that the parallax uncertainty is 
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largely dominating the final error on the linear diame- 
ter. Therefore, an intrinsically more precise interferomet- 
ric measurement would not result in an improvement of 
the linear diameter precision. The contribution of the 
VINCI/ VLTI measurement uncertainty in the error bar 
is only ±0.004 D©, or 37 % of the total error. 

5.4. Discussion 

Our LD angular diameter value can be compared to the 
value published by Hanbury Brown et al. ( 107 la). They 
found, using the Narrabri intensity interferometer, #ld = 
5.89 ± 0.16 mas. The difference with our value of 6*ld = 
6.039 ± 0.019 mas is -0.8 a, making the two results statis- 
tically compatible within their error bars. 

Davis & Tango (.1986) have obtained a value of 6>ud = 
5.63 ± 0.08 mas using the amplitude interferometer, and 
revised the Hanbury Brown et al. I|1974ajl value to 6*ud = 
5.60 ± 0.07 mas. Using a limb darkening coefficient u — 
0.5900 (B band from Claret et al. 2000) and the conversion 
factor formula from Hanbury Brown et al. I|1974b|l : 



Old = / l-u/3 
6ud ~ V 1 - 7«/15 



(8) 



we find p = 1.053. This translates into #ld = 5.93 ± 
0.08 mas and #ld = 5.90 ± 0.07 mas for the above UD 
values, respectively. They differ respectively by -1.3 and 
-1.9(7 from our result. This difference, marginally signifi- 
cant, could be explained by an uncertainty on the limb 
darkening factor. It is much stronger at visible wave- 
lengths than in the infrared, and can be affected by the 
presence of spectral features (Tango & Davis 2002 ). 

Using spectro-photometric observations, Cohen et 
al. ( 1992) have derived the angular diameter of Sirius and 
found 6.04 mas (with Kurucz model atmospheres), in re- 
markable agreement with our direct measurement. 

Due to the relatively slow rotational velocity of Sirius 
(■usini w 16 km.s -1 , from Royer et al. 12002ft. we do not 
expect any detectable flattening of its disk. In addition, 
the flux of Sirius B is totally negligible compared to A, 
in particular in the K band. We therefore do not foresee 
an asymmetry of the visibility function in azimuth due to 
these two contributors. 

Hanbury Brown et al. I|1974a[) have observed a small 
asymmetry of the Sirius visibility function, proposing ten- 
tatively that this may come from a dust disk around 
this star. In the error bars, our measurements do not 
show this asymmetry, as the D1-C3 and E0-G1 points are 
in good agreement with our single disk visibility model. 
Nevertheless, it should be noted that we have probed only 
two azimuth values, namely 80 and 160 degrees (D1-C3 
and E0-G1 baselines, respectively), and that our error bars 
on the diameter measurements on the shorter baseline are 
relatively large. Therefore, we cannot exclude formally the 
presence of a disk accounting for approximately 1% of the 
stellar flux. 



As a remark, Benest & Duvent 119951 have proposed 
that a late M5 dwarf could be orbiting Sirius A, but our 
measurements do not have the necessary sensitivity to 
record the visibility modulation that would be produced 
by a star of this type, with mn > 5. 

6. Modeling 

6.1. Hypothesis 

We computed a number of evolutionary models of Sirius A 
using the CESAM code (Morel \iWf^ including the pre- 
main sequence evolution. The ordinary assumptions of 
stellar modeling are made, i.e. spherical symmetry, no ro- 
tation, no magnetic field, no mass loss. The relevant nu- 
clear reaction rates are taken from the NACRE compila- 
tion (Angulo et al. EES) ■ The equation of state adopted 
is EFF (Eggleton et al. H973|l . and the OPAL opacities are 
from Iglesias & Rogers (1996) with the Grevesse & Noels 
(1993) mixture. The microscopic diffusion is described us- 
ing the formalism of Burgers ( 196SJ) with the resistance 
coefficients of Paquette et al. I|1986ll . We take into ac- 
count the radiative diffusivity as recommended by Morel 
& Thevenin (2002J), that limits the efficiency of the mi- 
croscopic diffusion in outerlayers of stars with intermedi- 
ate masses. We have neglected the radiative accelerations, 
as well as the changes of abundance ratios between met- 
als within Z. The atmosphere is restored using Hopf's 
T(t) law fMihalas H978|) . The adopted radius of the star 
is the bolometric radius, where r(r^) = T e ff. In the con- 
vection zones the temperature gradient is computed ac- 
cording to the MLTcm convection theory with a mixing 
length parameter of A = 1 (Canuto & Mazzitelli 119911 
1992). Following the prescriptions of Schaller et al. |1992) 
we have computed models that include overshooting of 
the convective core (radius R co ) over the distance O v = 
Amm(H p , R co ). We adopted a value A=0.15 in agreement 
with the results by Ribas, Jordi & Carme (2000) for a 
2 Mq star. In addition, we have also computed stellar mod- 
els without overshooting (A=0.00) of the convective core 
and with A=0.20 to estimate the effect of this parameter 
on the age of the star. As a comparison, we have also built 
a model with enhanced Z s /X, corresponding to the atmo- 
spheric abundances published by Qiu et al. I|2001J) . It is 
intended to demonstrate the strong effect of this overes- 
timated value of the internal abundance on the evolution 
stage of Sirius. Each model is described by about 760 mass 
shells. We considered that our model was representative 
of Sirius A when we reached the intcrferometric radius 
and the observed luminosity and effective temperature of 
this star (within their respective error bars). The main 
characteristics of the computed models are presented in 
Table E 

6.2. Results and discussion 

The HR evolution diagram is presented on Fig. 0] We 
adopt the standard definition of the ZAMS, i. e. where 
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Fig. 4. Evolutionary tracks in the HR diagram of models 
of Sirius A. The dashed rectangle delimits the uncertainty 
domain in luminosity and effective temperature, while the 
shaded area represents the uncertainty on the interfero- 
metric radius. The continuous line corresponds to model b 
with overshoot and the dashed-dot line to model a with- 
out overshoot, both with a mass of 2.12 M Q . Model d, 
corresponding to M=2.07 M Q and an age of 243 Myr, is 
represented by a dashed line. The filled square is the cen- 
ter of the rectangle. 



nuclear reactions begin to dominate gravitation as the pri- 
mary source of energy by at least 50%. For the mass of 
Sirius A this corresponds to an elapsed time of 4.2 Myr. 
Note that 100% of energy comes from nuclear reactions 
only after 10 Myr. The evolution of models starts with 
the quasi-static contraction of a cloud with a central 
temperature of 0.5 MK. According to these definitions, 
our model b of Sirius, which fits the VLTI/ VINCI radius 
and satisfies the adopted luminosity and T c g, is aged of 
197 ± 12 Myr. The error bar is fixed by the uncertainty on 
the VINCI/VLTI radius, that is itself largely dominated 
by the Hipparcos parallax error bar. 

Because the age of the WD can reach 210 Myr 
(Holberg, private communication) we also computed mod- 
els with a lower mass 2.07 MQ which fit both the observed 
radius and the T e g. We had to change sligthly Yi and the 
value of (y) [ within its error bar, see model d, resulting 
in an age of 243 Myr for Sirius. We can produce an derive 
a greater age of 250 Myr (model e) when increasing the 
overshooting parameter A to the value A=0.20. At the 
end of the pre-main sequence, Sirius A has generated a 
convective core due to the equilibration of the CNO nu- 
clear cycle. At the expected age of Sirius A for model d 
the hydrogen mass fraction at center is X c = 0.590 with 
a temperature of 22.46 10 6 K. The base of the external 
convective zone is situated at 0.9926 i?* and Y s = 0.200. 

Removing the overshooting from our model reduces 
the age of the binary by 10% (model a) while the model 
c with enhanced abundances leads to a too young age (25 
Myr) for Sirius A, in significant disagreement with the age 
of Sirius B. At the present status of our observational con- 



Table 3. Characteristics of our a CMa A models remain- 
ing within the uncertainty box of the observed properties 
of this star. The "Enhanced" column refers to the en- 
hanced abundances model. The subscripts i and s refer 
respectively to the initial content at zero age and the sur- 
face content at the end of the computed evolution. The 
adopted luminosity is log(L/L0) = 1.405 ± 0.022 with 
T e ff = 9900 ± 200K. The observed abundance of oxygen is 
-0.32±0.20 dex (see text). The precision of the age is ±12 
Myr corresponding to the crossing time of the shaded par- 
allelogram. 



Model 


a 


b 


c 


d 


e 


M/M Q 


2.12 


2.12 


2.12 


2.07 


2.07 


age (Myr) 


180 


197 


25 


243 


250 


overshoot 


0.00 


0.15 


0.15 


0.15 


0.20 


Toff (K) 


9 870 


9 906 


9 817 


9 936 


9 944 


log(L/L ) 


1.398 


1.403 


1.385 


1.408 


1.409 


log g 


4.30 


4.29 


4.30 


4.29 


4.29 


Yi 


0.267 


0.267 


0.330 


0.273 


0.273 




0.0165 


0.0165 


0.0350 


0.0145 


0.0145 


(*). 


0.0116 


0.0115 


0.0289 


0.0099 


0.0099 


[0/H] ; 


-0.17 


-0.17 




-0.23 


-0.23 


[0/H] s 


-0.33 


-0.33 




-0.39 


-0.39 


R/Rq 


1.7111 


1.7107 


1.7060 


1.7105 


1.7103 


Avq (mHz) 


82 


82 


82 


81 


81 



straints on Sirius A, it is impossible to decide if overshoot 
exists or not in the central convective core. Note that our 
model gives a surface abundance of [M/H]!, 11 ' 1 " 3 = —0.32 
dex which agrees with the abundance of the oxygen found 
by Kamp, Hempel & Holweger (2002). Finally, if Sirius A 
is able to excite radial modes of oscillation, we predict the 
mean large frequency spacing Ai>o to be ranging between 
81 and 82 fjHz This is the primary observable of aster- 
oseismology, and corresponds to the difference between 
the frequencies of oscillation modes with consecutive ra- 
dial order n (see Table |3J). This large spacing is sensitive 
primarily to the stellar density. It could therefore give a 
direct, high accuracy estimation of the mass of Sirius A 
when combined with the observed radius reported in this 
paper. 

7. Conclusion and perspectives 

We have reported in this paper our long-baseline interfer- 
omctric observations of a CMa A using the VINCI/VLTI 
instrument. In conjunction with constraints from previ- 
ous spectro-photometric observations, we derived an evo- 
lutionary model for Sirius A. We propose that the appar- 
ently high surface metallic content of Sirius is not charac- 
teristic of the whole average value of Z for the star, and 
is caused by the levitation of the heavy elements on the 
thin upper convective layer of Sirius A. By means of our 
model, we derive an age of 200 to 250 ± 12 Myr, consis- 
tent with the two evolutionary time estimates of Sirius B. 
The oldest age of Sirius A corresponds to a lower mass of 
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2.07 Mq which is the lowest mass value acceptable within 
its error bar. 

The accuracy on this age is greatly strengthened by 
the VLTI/ VINCI radius, thus encouraging further stud- 
ies to improve our knowledge of the diameter of nearby 
stars. Based on our model, we predict that the asteroseis- 
mic large frequency spacing of Sirius A should be 81 to 
82 /iHz, if it exhibits radial oscillations. A more complete 
modeling of the external layers of this star for asterosismic 
frequencies prediction, including radiative accelerations, is 
currently in preparation. 

We would like to emphasize here the need for very high 
accuracy parallax values, beyond the Hipparcos preci- 
sion. We have shown that the uncertainty on our determi- 
nation of the linear radius of Sirius A is largely dominated 
by the error on the distance, despite its proximity (2.7 pc). 
The next generation of interferometric instruments, such 
as the AMBER beam combiner (Petrov et al. 2000J), will 
provide high precision angular diameters, together with 
their wavelength dependence. In order to constrain the 
stellar structure models, that output linear values, it will 
be necessary to match the precision on the distance to the 
angular diameter uncertainty. It is expected that future 
space based astrometric missions will be able to give a 
significant improvement over the Hipparcos parallaxes, 
and are therefore crucial in this respect. This work on 
Sirius A and the study of a Cen A & B by Kervella et 
al. l|2003bjl have demonstrated that for the nearest stars, 
thanks to the very small error bar on the linear radius, we 
have advantage to replace the classical error bars in the 
HR diagram, e.g. luminosity and effective temperature, by 
the couple luminosity and radius. 

The direct measurement of the limb darkening of 
Sirius A is the next step of the interferometric study of 
this star. It will allow a refined modeling of its atmosphere 
and will complete the calibration of the radius that we 
presented in this paper. This observation will be achieved 
soon using the long baselines of the VLTI. 
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